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Calculations of Potential Energy Surfaces Using Monte Carlo Configuration 
Interaction 

Jeremy P. Coe,-"^ Daniel J. Taylor/ and Martin J. Paterson^'F^ 

Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, 
EHI4 4AS, UK. 

We apply the method of Monte Carlo configuration interaction (MCCI) to calculate ground-state potential 
energy curves for a range of small molecules and compare the results with full configuration interaction. We 
show that the MCCI potential energy curve can be calculated to relatively good accuracy, as quantified using 
the non-parallelity error, using only a very small fraction of the FCI space. In most cases the potential curve 
is of better accuracy than its constituent single-point energies. We finally test the MCCI program on systems 
with basis sets beyond full configuration interaction: a lattice of fifty hydrogen atoms and ethylene. The 
results for ethylene agree fairly well with other computational work while for the lattice of fifty hydrogens we 
find that the fraction of the full configuration interaction space we were able to consider appears to be too 
small as, although some qualitative features are reproduced, the potential curve is less accurate. 
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I. INTRODUCTION 

Full configuration interaction (FCI) enables, in theory, 
a quantum system to be modelled as accurately as pos- 
sible within a given basis set, but the rapidly increasing 
number of configurations means that such calculations 
are out of reach except for sufficiently small molecules 
and basis sets. Standard truncation methods, such as 
considering only single and double excitations with re- 
spect to a reference state or restricting excitations to a 
pre-defined space, can reduce the size of the calculation 
but often at the expense of the accuracy and consistency 
of the correlation energy as important configurations may 
be neglected. 

It is acknowledged, however, that a large proportion 
of the states in many full configuration interaction wave- 
functions tend to have practically negligible coefficients 
and should not be expected to individually contribute 
substantially to the properties of the wavefunction. Novel 
approaches to truncated configuration interaction have 
been developed to attempt to tackle the difficulties posed 
by the exceedingly rapid growth of the standard full CI 
space by exploiting this observation and trying to seek 
out the important states, see Ref. fl. for a recent review. 
For example these include a priori estimates^i^ Monte 
Carlo sampling of a DMRG calculation'* and estimates 
from perturbation theory.^ Another promising method is 
that of Monte Carlo configuration interaction (MCCI) 
developed originally by Greer<^i In MCCI a trial wave- 
function is randomly augmented with coupled states in 
an iterative scheme where those states with coefficients 
smaller than a certain value in the resulting solution of 
the Schrodinger equation are eventually removed. Even 
without prior knowledge of the important configurations 
or molecular orbitals such a procedure can, in principle, 
result in a compact wavefunction which recreates much 
of the energy of the FCI wavefunction but using only a 
small fraction of the FCI states. 
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Single- reference methods based on coupled-cluster— 
such as coupled-cluster singles doubles (CCSD)2. are con- 
sidered to have some of the best efficiencies in calcula- 
tions for systems in which the FCI expansion would have 
one dominant configuration and the correlation is con- 
sidered dynamic. However, they may have problems — 
especially if the restricted Hartree-Fock determinant is 
used as a reference — when there are a number of impor- 
tant configurations and the system is deemed multirefer- 
ence. These configurations are associated with static cor- 
relation and, e.g., dissociation energies may be modelled 
incorrectly. For multireference systems the associated 
static correlation may be accounted for by using meth- 
ods such as complete active space self-consistent field 
(CASSCF)ji2, but insight is required into the selection of 
orbitals for the active space and the calculation becomes 
intractable as the size of the active space increases. If 
the orbitals have been selected appropriately then per- 
turbation methods (CASPT2) or multi-reference config- 
uration interaction (MRCI) may then be used to account 
for the remaining dynamic correlation but at relatively 
high cost. In principle MCCI can give a compact wave- 
function which can account for both static and dynamic 
correlation with minimal user choices: the accuracy and 
size of the calculation is essentially controlled only by the 
coefficient cut-off parameter Cmin- 

MCCI has been successfully applied to the single point 
energy of the water molecule)^ and the dissociation en- 
ergy of HF and H20.*"'^ Electronic excitation energies 
have also been calculated*^ for the first-row atoms begin- 
ning with carbon — and also for silicon — with relatively 
small average errors and numbers of states. While the 
excitation energies of molecules CII2, C2, N2 and H2O 
were considered in Ref. [l^ with only tens of meV error 
and using only a very small proportion of the states com- 
pared with FCI. 

In this work we investigate the ability of MCCI to cal- 
culate potential energy curves. If the MCCI results can 
have an almost balanced error across different geometries 
then we should be able to produce a curve of better ac- 
curacy, as quantified by the non-parallelity error, than 
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its constituent single-point energies. We first assess the 
usefulness of this method on molecules with basis sets for 
which FCI results exist. We look at potential curves with 
respect to the dissociation of one hydrogen from HF, BH, 
and CH4. The curves for bond stretches are produced for 
C2, F2, N2 and H12. The potential curves for NH3 inver- 
sion and the model formation of BeH2 along the reaction 
coordinate are also calculated. Finally systems beyond 
the scope of current FCI are also investigated. We cal- 
culate a potential curve for H50 which we compare with 
density-matrix renormalization group (DMRG)i^ results, 
and a potential curve for ethylene as the torsional angle 
is varied which we compare with other computational re- 
sults. 



II. METHOD 

A FCI wavefunction may be written in the notation of 
second quantization as 

|vI')=co|*o) + E^«]«''l*o) 

+ J2 4«Jafea]«d*o) + ■ • • (1) 
k<i.l<j 

where (a^) creates (annihilates) the orbital i in a 
state. Here letters appearing as subscripts on the coefB- 
cients label occupied orbitals in the reference state |^o} 
while superscripts are unoccupied orbitals in the refer- 
ence. A traditional truncation would limit the calcula- 
tion size by only including substitutions up to a certain 
substitution level. Limiting the expansion to only the de- 
picted terms in EqUwould correspond to a CISD calcula- 
tion. MCCI instead randomly augments a wavefunction 
4'} = '^^Ci lipi) with coupled configurations, i.e., states 
{ipj I such that {ipj \H l'^) ^ so that, in principle, the im- 
portant configurations can eventually be found regardless 
of their substitution level. Configuration state functions 
(CSFs) are used in MCCI rather than the more usual 
Slater determinants (SDs). This means that the MCCI 
wavefunction fulfils at least one of the requirements of the 
exact wavefunction: it is an eigenfunction of 5*^. In addi- 
tion fewer states are needed, e.g., the ratio of SDs to the 
number of CSFai^ when S = Mg = for 8 electrons and 
30 basis functions is around 4.4. One downside is that 
the calculation of the overlap and Hamiltonian matrix is 
more complicated than when using SDsJ^ 

We briefiy describe the full procedure for the MCCI 
algorithm^'^ which consists of branching, diagonalization 
and pruning. In the branching step, coupled CSFs, cre- 
ated by a symmetry preserving stochastic single or dou- 
ble substitution with equal probability in the current set 
of CSFs, are added. Branching is always attempted from 
CSFs with a coefficient greater than a certain value while 
other CSFs each have probability one half of being used. 
The coefficients are then found for the approximate wave- 
function expanded in this augmented set of CSFs by solv- 



ing the generalised eigenvalue problem of He = ESc. 
Here H and S are the Hamiltonian and overlap matrix 
within this subspace. Finally there is a pruning step: 
added CSFs with coefficients less than a certain value, 
Cniin, in the wavefunction expansion are then discarded. 
These three steps are repeated for a large enough num- 
ber of iterations so that the energy appears to have con- 
verged. In addition, every k iterations (where fc = 10 
in this work) a pruning step is implemented whereby all 
of the CSFs in the wavefunction are also considered as 
candidates for removal depending on the magnitude of 
their coefficient. Furthermore there is no branching or 
pruning on the last iteration, but a full pruning step is 
carried out on the penultimate iteration. The program 
can run in parallel where states that have been added 
and retained at each step on each processor are shared 
with all other processors. The results in this work come 
from MCCI calculations using either eight or twelve pro- 
cessors. 

We initialise the procedure with the CSF formed from 
the restricted Hartree-Fock wavefunction within a given 
basis. New Hartree-Fock orbitals are calculated at each 
geometry. For the calculation of the Hartree-Fock molec- 
ular orbitals and their one-electron and two-electron in- 
tegrals we use the MOLPRO packagcjii 

III. RESULTS 

As the potential is defined up to an additive constant 
then it is the shape of the potential energy curve that is 
important. So two curves that differed only by a constant 
would for all practical purposes be the same. A useful 
approach to quantify the accuracy of the MCCI poten- 
tial curve is therefore the non-parallelity error (NPE)i^ 
defined here for an approximate energy E^^'p™^ as 

NPE = max I^J^d _i;|^^P'°-| _nun -£;Jp™'^| (2) 

where R ranges over all considered bond lengths. 

Hence if we can achieve a balanced accuracy in the 
energy across the geometries we can recover a potential 
energy curve with higher accuracy than its single-point 
energies. We note that this is not necessarily guaranteed 
as if many points are essentially the FCI energy but there 
are a few points with large errors then the NPE would be 
high but the mean single-point energy error could be rel- 
atively low. To investigate this we attempt to use a small 
enough Cmin, where computationally feasible, to produce 
a sufficiently accurate curve. 

A. Hydrogen dissociation 

We first consider varying the length of a bond with hy- 
drogen in BH, HF, and CH4 and compare the potential 
energy curves with FCI results from Ref. 19'. These po- 
tential curves represent one of the simplest dissociations: 
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the breaking of a bond to hydrogen. As such they would 
be expected to have some multireference character and 
make an interesting test case for MCCI to be compared 
with other methods. For these three systems one core 
orbital is frozen in all the results. 



Figures [2] and |3] show that the MCCI results are clos- 
est to the FCI when the bond length is large. We also 
see that rerunning the calculations with a different ran- 
dom number seed does not change the error much and 
the variation in error seems to be lower as the cut-off is 
decreased. 



1. HF 

We see in Fig. [T] that the shape of the HF FCI curve 
in a 6-31G* basis is reproduced by MCCI and it appears 
that for this system, apart from the energy shift, there is 
not much difference between the two cut-off values. The 
reduction in NPE by decreasing the cut-off value from 
Cmin = 5 X 10~^ to Cmin = 5 X 10"'* is 5.7 to 1.3 kcal/mol. 
These NPE compare very favourably with methods based 
on a single reference where CCSD gives a value of 13.0 
kcal/mol, the CCSD(T) NPE is as high as 26.8 kcal/mol, 
while that of UCCSD(T) is 3.7 kcal/mol.i^ For the mul- 
tireference approaches considered in Ref. i20l, an NPE of 
17.69 kcal/mol was found for CASSCF(8, 5) 'valence ac- 
tive space' and 4.83 kcal/mol for CASSCF(8, 8) '1:1 ac- 
tive space'. CASPT2(8, 5) and CASPT2(8, 8) resulted in 
NPEs of 2.77 and 0.5 kcal/mol. The NPEs for S0CI(8, 5) 
and S0CI(8, 8) were 3.20 and 0.04 kcal/mol respectively. 
Ref. |21| looked into using a minimal active space (one 
bonding orbital and one anti-bonding orbital for each 
broken bond) and although their CASSCF(2, 2) had a 
higher NPE (18.66 kcal/mol) than when using the afore- 
mentioned active spaces, the CASPT2 value was lower at 
0.47 kcal/mol and the SOCI was better than when using 
the larger valence active space but, at 0.98 kcal/mol, not 
as good as when using the 1:1 active space. 

The larger cut-off for MCCI required 172.3 CSFs on 
average, while the smaller needed 1,337. The full con- 
figuration space consisted of 3, 756, 816 SDs^^ so only a 
very small fraction were used by MCCI to give an NPE 
of a few kcal/mol (see table IJ). 
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FIG. 1. MCCI and FCI results^- for energy (Hartree) against 
bond length R (angstrom) for HF with the 6-31G** basis set 
and one frozen core orbital. 
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FIG. 2. Error in two MCCI runs with Cmin = 5 x 10"^ 
compared with the FCI results^^ against bond length R 
(angstrom) for HF with the 6-31G** basis set and one frozen 
core orbital. 
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FIG. 3. Error in two MCCI runs Cmin = 5 x 10"'' compared 
with the FCI results^'' against bond length R (angstrom) for 
HF with the 6-31G** basis set and one frozen core orbital. 



As the MCCI potential curves (Fig. [TJ appear to have 
essentially reached dissociation, we may compare MCCI 
energies at i? = 4 angstrom with the total energy from 
appropriate MCCI calculations on the fragments to in- 
dicate how close the MCCI calculation is to size consis- 
tency. We find a difference of 4.0 kcal/mol for the larger 
cut-off and 0.15 kcal/mol for the smaller cut-off. 
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2. BH 

The results for BH were calculated in an aug-cc-pVQZ 
basis. Preliminary MCCI results for BH are contained 
in Ref. (l). In Fig. U we can see that with a cut-off 
of 5 X 10"'^ the general shape of the MCCI potential 
curve is qualitatively correct with the minimum in ap- 
proximately the right place and a well-behaved dissocia- 
tion although there is a small shoulder in the potential 
around R — 3 angstrom. The limiting value of the en- 
ergy is reached too early however. By decreasing Cmin 
to 5 X lO-'^ the curve is almost that of the FCI shifted 
by a small energy. The NPE is 22.8 kcal/mol in the for- 
mer case and reduces to 2.6 kcal/mol with the smaller 
Cmin value. For this value of Cmin the mean number of 
CSFs was only 4, 220 compared with around fifteen mil- 
lion Slater determinants in the FCI space (See Table |l| 
while 330 were required on average for the larger case. 
From Ref. 1^ the NPE value for CCSD was 8.1 kcal/mol, 
that of CCSD(T) was 23.3 kcal/mol while UCCSD(T) 
had an NPE as low as 3.1 kcal/mol. Multireference 
methods have also been considered^^ and compared with 
FCI results for the potential curve of BH but with a 
smaller range and a cc-pVQZ basis as there were diffi- 
culties in convergence when using the aug-cc-pVQZ ba- 
sis. Although not a direct comparison, we note that the 
CASSCF(4, 4) 'valence active space' there had an NPE 
of 12.68 kcal/mol, CASSCF(4, 5) '1:1 active space' gave 
9.38 kcal/mol, CASPT2(4, 5) gave 3.16 kcal/mol while 
only the second-order CI (SOCI) results were lower than 
the MCCI value in this work at 0.29 and 1.54 kcal/mol 
for active spaces of (4, 4) and (4, 5) respectively. 
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FIG. 4. MCCI and FCI results^^ for energy (Hartree) against 
bond length R (angstrom) for BH with the aug-cc-pVQZ basis 
set and one frozen core orbital. Adapted from Ref. [l|. 

We see in Fig. [5] that the energy error for this sys- 
tem at Cmin = 5 X 10~^ when compared with FCI gen- 
erally decreases as the bond length increases. We note 
that the difference between the MCCI energy at i? = 6 
angstrom and the total MCCI energy for the fragments 
is 6.6 kcal/mol at Cmin = 5 x 10~^ and 0.83 kcal/mol at 



Cmin — 5 X 10 ^. 




R (angstrom) 

FIG. 5. MCCI (cmin = 5 X 10"*) energy error (Hartree) 
when compared with FCI results against the bond length R 
(angstrom) for BH with the aug-cc-pVQZ basis set and one 
frozen core orbital. The NPE is 2.6 kcal/mol. 



3. CH^ 

As in Ref. fl9' a tetrahedral geometry was used for 
methane with bond length 1.086 angstrom for the three 
carbon-hydrogen bonds that were not varied. Prelim- 
inary MCCI results for CH4 are contained in Ref. [l|. 
Fig. [6] shows that, in a 6-31G* basis, MCCI again cap- 
tures the shape of the FCI potential curve. Now it ap- 
pears that the larger cut-off may be a little too high in 
energy at large R compared with the equilibrium geom- 
etry. This is quantified using the NPE which is 10.3 
kcal/mol for Cmin = 5 x 10"'^ but drops to 0.6 kcal/mol 
at the smaller cut-off. We note that in Ref. [H CCSD 
gave a value of 10.3 kcal/mol, UCCSD gave 5.1 kcal/mol 
while UCCSD(T) was as low as 3.2 kcal/mol. In Ref. [13 
the definition of the 1:1 active space was equivalent to the 
valence active space in this case and gave CASSCF(8, 8) 
with an NPE of 6.34 kcal/mol, CASPT2(8, 8) gave 1.56 
kcal/mol, CISD[TQ](8, 8) gave 1.33 kcal/mol, while only 
S0CI(8, 8) had a lower NPE than the MCCI resuhs with 
0.3 kcal/mol. The minimal active space results^! had an 
NPE of 9.25 kcal/mol for CASSCF(2,2), with CASPT2 at 
1.22 kcal/mol and SOCI a little worse with 0.6 kcal/mol. 
Hence only one of the results from these two works using 
multireference methods gave a lower NPE for methane. 

The FCI space is around 26.7 million SDs,^^ but in 
Table U we see that the MCCI results are achieved using 
only 417 and 4, 272 CSFs on average for the larger and 
smaller Cmin respectively. 

In Fig. [7] we see that the error in the energy compared 
with FCI is greatest at very small bond lengths and now 
the minimum error is at bond lengths a little longer than 
that at equilibrium, but the errors are all very small. 
The difference in energies between MCCI applied to the 
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FIG. 6. MCCI and FCI resultsi^ for energy (Hartree) against 
one carbon-hydrogen bond length R (angstrom) for CH4 with 
the 6-31G* basis set and one frozen core orbital. Adapted 
from Ref . Q. 



TABLE I. MCCI mean CSFs compared with FCI symmetry- 
adapted SDs 



Method 



BH states HF states CH4 states 



MCCI0.005 (CSFs) 
MCCI0.0005 (CSFs) 
FCI (SDs) 
MCCI0.005 fraction 
MCCI0.0005 fraction 



333.1 
4219.6 
15,132,412 
0.0022% 
0.028% 



172.3 
1337.0 
3, 756, 816 
0.0046% 
0.036% 



417.0 
4272.0 
26, 755, 625 
0.0016% 
0.016% 



than many of the multireference results from the litera- 
ture. Some CASPT2 and SOCI calculations, using large 
active spaces or chemical intuition in the construction of 
a smaller active space, produced a lower non-parallelity 
error but the number of states used for these calculations 
was not reported for comparison. 



system at i? = 4.6 angstroms and the sum of MCCI ap- 
plied to the fragments is 15 kcal/mol at the larger cut- 
off and 0.79 kcal/mol at the smaller cut-off. For HF, 
BH and CH4 at Cmin = 5 x 10~^ the difference between 
the FCI and MCCI energy, at all geometries considered, 
ranged from around 8.4 kcal/mol to 33.8 kcal/mol while 
this range was around 1.1 kcal/mol to 4.5 kcal/mol when 
the cutoff was lowered to Cmin = 5 x lO"**. 
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FIG. 7. MCCI (cmin = 5 X 10"") energy error (Hartree) when 
compared with FCI results against one carbon-hydrogen bond 
length R (angstrom) for CH4 with the 6-31G* basis set and 
one frozen core orbital. The NPE is 0.6 kcal/mol. 

We have seen, for the three systems considered for hy- 
drogen dissociation, that by going to a low enough co- 
efficient cut-off value of Cmin = 5 X 10~^ we can repro- 
duce the shape of the potential curve to relatively high 
accuracy: we achieve a non-parallelity error of around 
a few kcal/mol which is better than the single- reference 
methods we compare with from the literature. This is 
accomplished using in the order of hundredths of a per- 
cent of the full CI space. The accuracy was also higher 



B. Carbon Dimer 

Next we look at the dissociation of C2 in a 6-31G* basis 
with two frozen orbitals. This system is known to be 
multireference and possess low-lying excited states hence 
poses a stern test of MCCI. Preliminary MCCI results 
for C2 are contained in Ref. [l|. FCI results^S and the 
MCCI potential curve with Cmin = 5 x 10~^ are shown in 
Fig. [5] The MCCI curve accuracy is somewhat difficult to 
judge by eye due to the many close excited states, but it 
appears to follow the FCI ground state closely except for 
bond lengths (R) around 3 angstrom where it seems that 
a low lying excited state (B') has been converged upon 
instead. However it appears that the three FCI curves 
are in fact tending to a degenerate state as the system 
dissociates. In addition, the lowest energy state between 
about 1.7 and 2.5 angstrom is B rather than X and MCCI 
may have converged to this as would be expected: the 
states have different symmetry so can cross but have the 
same spatial symmetry in the abelian subgroup {D2h) 
used for computation. 

The mean number of CSFs employed was ~ 6, 900 
for MCCI here, while FCI results from Ref. [23 re- 
quired 52,407,353 symmetry-adapted Slater determi- 
nants. This gives a ratio of 0.013%. The NPE is 
4.9 kcal/mol for MCCI here which is substantially bet- 
ter than methods based on a single reference. From 
Ref. m CCSD had an NPE of 24.3 kcal/mol and that 
of CCSD(T) was 61.3 kcal/mol, highlighting the insta- 
bility of coupled cluster at large R. Unrestricted meth- 
ods behaved better at dissociation but were less good 
at intermediate bond lengths resulting in a NPE for 
UCCSD(T) of 21.6 kcal/mol. CISDTQ gave a curve 
with 16.6 kcal/mol for the NPE. It is worth noting that 
a CISDTQ, without symmetry, and using Slater deter- 
minants would be expected to have around three mil- 
lion terms. Even allowing for the reduction in the size 
of the space through spatial symmetry and the use of 
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FIG. 8. MCCI (cmin = 5 X 10"'') and FC& results for energy 
(Hartree) against the bond length R (angstrom) for C2 using 
the 6-31G* basis set with two frozen core orbitals. Here X 
and B' are ^Ei" and B is ^Ag. Adapted from Ref. [ll. 



CSFs, the MCCI result with circa 6, 900 CSFs on aver- 
age and a better NPE shows the usefulness of the MCCI 
approach here. Valence active space CASSCF in Ref. 
used 660 determinants and had an NPE of 5.4 kcal/mol 
for (SA)-CASSCF. (EOM)CCSD gave 24.3 kcal/mol and 
CR-(EOM)CCSD(T),III gave 13.5 kcal/mol for the NPE. 
CISD[TQ] had an NPE of 16 kcal/mol in that work 
and used 87, 415 determinants. In fact only two of the 
multireference techniques investigated in Ref. [2^ had a 
lower NPE: (SA)-CASSCF MRCI with 0.4 kcal/mol us- 
ing 270, 338 determinants and SA-CASPT2 with an NPE 
of 3.8 kcal/mol. 

A wavefunction for the carbon dimer found using FCI 
quantum Monte Carlo (FCIQMC) in Slater determi- 
nant space in Ref. 24 required 9.15 x 10^ walkers at 0.9 
angstrom in this 6-31G* basis to give almost FCI results 
although the point at i? = 3 angstrom was neglected. We 
note that with a cc-pVDZ basis and requiring that only 
states with zero angular momentum were allowed they 
calculated very similar results to MRCI, and with lower 
energy for R > 1.15 angstrom than when using 6-31G*, 
with 2 X 10^ walkers compared with a standard FCI space 
of 4.74 X 10^ Slater determinants. This suggests that the 
performance of MCCI here could perhaps be improved if 
we also discriminate between states using their angular 
momentum in addition to their symmetry and spin. 

We note that MCCI potential curve for the carbon 
dimer is of slightly higher accuracy, in a sense, than its 
constituent single-point energy calculations as the mean 
single-point error was around 6.0 kcal/mol. In addition, 
if we exclude the two points around i? = 3, that seem to 
have converged to a low-lying excited state, then the sin- 
gle point energy error only reduces a little to 5.7 kcal/mol 
but the NPE is now just 2.7 kcal/mol. In Fig. [9] we see 
that the error does not vary too much, as summarised 
by the NPE, but is highest approaching the largest bond 
length. Interestingly it is lowest not at equilibrium but 



at bond lengths from around 1.8 to 2 angstroms, which 
in this case may be due to the lowest energy state be- 
ing reached, not the original ground-state. This slightly 
lowers the mean single-point error in this case, and will 
slightly raise the NPE. At 3 angstrom we compare the 
MCCI energy for the full system with that of the sum of 
the MCCI energies of the fragments and find a difference 
of 3.2 kcal/mol. We note that as the energy of the system 
using MCCI is higher than the energy of the fragments 
here then as the system has perhaps not sufficiently dis- 
sociated the difference is perhaps smaller than if a larger 
R had been considered. 



CO 

X 
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FIG. 9. MCCI (cmin = 5 X 10"*) energy error (Hartree) 
when compared with FCl2^ results against the bond length R 
(angstrom) for C2 using the 6-31G* basis set with two frozen 
core orbitals. The NPE is 4.9 kcal/mol. 



1. Analysis of the MCCI wavefunction 

At the equilibrium geometry for C2 the MCCI wave- 
function has states with as many as 7 substitutions from 
the Hartree- Fock reference and the multiconfigurational 
nature is evident from the largest ten coefficients of the 
CSFs in Table im By i? = 2.0 the dominant configuration 
is now a double substitution of the Hartree- Fock reference 
which could suggest that the lowest lying state here (B), 
not the equilibrium ground-state curve (X), may have 
been converged upon as expected or could just be an as- 
pect of the multiconfigurational nature which is is more 
apparent here; the largest two coefficients are —0.636 
and 0.557. By R — 3.0 the largest coefficient belongs 
to a CSF which differs to the major contributor in the 
previous two cases. This may be linked to the suggestion 
from the potential curves that here a different excited 
state is converged upon at large R. 
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TABLE II. Ten largest coefficients of the C2 MCCI wavefunc- 
tion at i? = 1.25 angstrom with the number of substitutions 
with respect to the Hartree-Fock reference hsted by spin. 



Coefficient a substitutions (5 substitutions 

-0.830 

0.331 1 1 

-0.184 1 1 

-0.180 1 1 

-0.179 1 1 

-0.160 2 

0.158 1 1 

0.158 1 1 

-0.142 1 

0.106 2 2 



C. F2 

In Fig.[Tn]we now display the MCCI potential curve for 
F2 in a cc-pVDZ basis set without any frozen orbitals. 
The MCCI results can be seen to follow the shape of 
the full valence CI results of Ref. Hi despite the MCCI 
calculation not freezing any orbitals. We can also see 
how the Hartree-Fock energy increases much too rapidly 
as the bond length becomes large. 
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FIG. 10. MCCI (two frozen orbitals) with Cmin = 5 x 
10~*, Full valence CI^ and Hartree-Fock results for energy 
(Hartree) against bond length R (angstrom) for F2 using the 
cc-pVDZ basis set. 

We note that the MCCI with Cmin = 5 x 10~^ required 
around 3, 577 CSFs on average for the depicted points. 
The FCI space with two frozen orbitals would be ex- 
pected to comprise 4.3 x 10^^. Even though only a very 
small fraction of this space is used (8.3 x 10"^), when 
compared with the Full valence CI points the NPE for 
MCCI is 6.2 kcal/mol. 

In Fig.[TT]we display the difference between the MCCI 
and Full valence CI for the latter's points in the range 
displayed in Fig. [TUl Here the MCCI error tends to de- 



crease with increasing bond length for the points dis- 
played. Comparing the MCCI energy of the system at 
i? = 3 angstrom with the sum of the MCCI energies of 
the constituent atoms we find a difference of 7.3 kcal/mol. 
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FIG. 11. Energy error of MCCI (two frozen orbitals) with 
Cmin = 5 X 10~* when compared with Full valence CI^^ against 
bond length R (angstrom) for F2 using the cc-pVDZ basis set. 
The NPE is 6.2 kcal/mol. 



D. N2 

The nitrogen molecule is well known for exhibiting 
multireference character at long bond lengths and here 
methods based on a single reference can fail, see, e.g., 
Ref. [ii. Hence this is another challenging system to 
test if MCCI can produce a potential energy curve with 
balanced error. The MCCI potential curve for a cc- 
pVDZ basis with two frozen core orbitals is displayed in 
Fig. HH and compared with FCI results from a number of 
sourcesi^"— We see that the general shape of the curve is 
recovered by MCCI when using a cut-off of Cmin = 10^'^ 
and the gap appears fairly constant except at very small 
bond lengths. For the values for which FCI results could 
be found, depicted in Fig. [T^l the NPE is 6.6 kcal/mol, 
while the MCCI CSFs ranged from about a thousand to 
almost five thousand as the bond length increased with a 
mean value of 2, 854. This is in contrast to a SD space of 
4.3 X 10^ for a FCI when ignoring spatial symmetries. The 
200 iterations required from around 3 minutes at small 
R to 1.3 hours as R became large using 8 processors. 

For the six bond lengths considered in Ref. [2^ for 
this system using the density matrix renormalization 
group (DMRG) the NPE in mEh ranged from 0.4 to 
0.01 for DMRGIOOO and DMRG4000 respectively while 
CCSD(T) was 172.724 and CCSD 57.754. Using a 
CAS(6, 6) wavefunction as a reference, MRCI and MRCC 
had 0.464 and 0.732 mEh respectively, but the number 
of states used was not specified. While the results of 
our MCCI calculation with Cmm = 10"'^ gave 5.37mEh 
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FIG. 12. MCCI (cmin = 10"^) and FCI^i"^ energy (Hartree) 
against bond length 7? (angstrom) for N2 using the cc-pVDZ 
basis set with two frozen cores. 



but using only a few thousand CSFs on average. The 
variational parameters of DMRG allow a rough compar- 
ison with the size of the space. The number of param- 
eters is known to be 0{M^k) (see for example Ref. Isol ) 
and so for the 26 orbitals non-frozen orbitals used there 
are of the order of 26 million parameters for DMRGIOOO 
(M = 1000). Ref. [31I uses the same set of points to 
compare the accuracy of multireference coupled cluster 
approaches, again with a CAS (6, 6) reference. For state- 
specific multireference CCSD (SS-MRCCSD) with delo- 
calised orbitals they find an NPE of 12.0 mEh, single- 
root MR Brillouin-Wigner CCSD (srMRBWCCSD) had 
an NPE of 19.4 mEh, and MR averaged quadratic CC 
MRAQCC gave 1.1 mEh- 

In Ref. [2l|, the thirteen points of Ref. ^2?' are used. Over 
these points we calculate the MCCI NPE as 5.9 kcal/mol. 
The CASSCF(10, 10) '1:1 active space' calculation^ had 
the highest NPE at 22.93 kcal/mol of the methods they 
considered for this system, while the CASSCF(10, 8) 'va- 
lence active space' result had an NPE of 15.03 kcal/mol. 
Using their minimal active space, defined with knowledge 
of which bonds were broken, resulted in a CASSCF(6, 
6) with a NPE of 14.59 kcal/mol. Interestingly the 
CASPT2(6, 6) had the largest NPE of the three active 
spaces considered for CASPT2 at 5.2 kcal/mol, unlike 
their results for HE and cH4, while CASPT2(10, 10) had 
the smallest at 1.88 kcal/mol. The SOCI were all be- 
low one kcal/mol but it was not clear how many states 
were used to achieve this result. MCCI with the cut- 
off employed is hence relatively close in accuracy to the 
CASPT2 results where knowledge of which bonds are 
broken has been used and it appears that large active 
spaces for CASPT2 are needed to reduce the NPE to 
below 2 kcal/mol. 

In Fig. [13] we see that the difference between the MCCI 
and FCI energy tends to increase with the bond length 
for the nitrogen molecule. 

We have also calculated a potential curve using a cc- 
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FIG. 13. Error (Hartree) in the MCCI (c^in = 10~^) cal- 
culation when compared with the FCI against bond length 
R (angstrom) for N2 using the cc-pVDZ basis set with two 
frozen cores. The NPE is 6.6 kcal/mol. 



pVTZ basis. This is displayed in Fig. [14] and shows 
how the dissociation is again well-behaved for this larger 
range of points. Without symmetry considerations the 
FCI space is around 10^^ SDs while the mean number of 
CSFs required for the MCCI curve was around 4, 600 at 
Cmin = 10~^ and around 10, 600 at c^in = 5 x 10~^ . The 
very large space precludes a comparative FCI calculation. 
However we may make an approximate comparison with 
FCI results^ which used an ANO[4s3pld] basis with two 
frozen cores. The difference in energy between the MCCI 
results at the equilibrium bond length of 1.098 angstrom 
and at 4 angstrom is 0.381 Hartree at the smallest cut-off 
used. This compares reasonably well with the approxi- 
mate dissociation energy FCI results of Ref. '32I. There 
an approximate equilibrium bond length of 2.1 Bohr was 
used to give 0.321 Hartree and 40 Bohr was used to cal- 
culate the energy at dissociation. When comparing the 
MCCI energy of the system at i? = 4 angstrom with 
the sum of the MCCI energies of the appropriate frag- 
ments we find a difference of 56.8 kcal/mol when using 
Cmin = 10^'^ and a difference of 41.8 kcal/mol when us- 
ing Cmin = 5 X 10~^. This suggests that the atoms are 
treated more accurately when using MCCI with HE MOs 
on this multireference system with a basis that is beyond 
current FCI. However we note that the curve appears to 
dissociate correctly even with the largest difference seen 
in this work between the molecule essentially at dissocia- 
tion and the fragments suggesting that MCCI is the least 
size-consistent in this case. 

To give a somewhat fairer comparison with the FCI 
results of Ref. 113 we remove the f orbitals from a cc- 
pVTZ basis, and attempt to use the same geometries as 
Ref. but keep the MCCI cut-off value the same and 
do not freeze any core orbitals. For R = 2.1 Bohr we 
find an energy of —109.296 Hartree using 12,664 CSFs, 
but the MCCI results were much too high in energy at 
R = 40 Bohr (-93.66 Hartree using 5, 029 CSFs) which 
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FIG. 14. MCCI energy (Hartree) against bond length R 
(angstrom) for N2 using the cc-pVTZ basis set with a cut-off 
of Cmin = 10~^ Or Cmin = 5 X 10~*. Inset: Energy difference 
between the MCCI results using the two cut-offs. 



we attribute to the restricted Hartree-Fock (RHF) ref- 
erence being qualitatively wrong. The amount of states 
needed to compensate for this and produce a sensible 
MCCI energy using the RHF molecular orbitals is then 
perhaps much too large for the cut-off employed. How- 
ever we saw that when using a full cc-pVTZ basis the po- 
tential curve dissociated as one would expect in Fig. [T3] 
when the largest bond length was 4 angstrom. Hence we 
note that the problems resulting in a very poor reference 
at extremely large bond length could be circumvented, 
to a degree, by calculating the curve until approximate 
convergence towards the energy as dissociation where the 
smaller bond distance means that RHF is not quite as in- 
correct. To this end we use R — 7.56 Bohr (4 angstrom) 
and find the energy to be -108.959 with 19,156 CSFs. 
This gives an approximate dissociation energy of 0.337 
Hartree. The MCCI result therefore has an approximate 
error of around 10 kcal/mol. The symmetry adapted 
FCI space was around 9.68 x 10^ while the mean number 
of CSFs used by MCCI for the two points was 15, 910. 
Hence here only around 1.6 x 10^'*% of the FCI symmetry 
adapted space was used. 



glecting spatial symmetries, while the MCCI wavefunc- 
tion with Cmin = 10~^ consisted of CSFs ranging from 386 
at a; = to a maximum of 1,002 around the transition 
state at X = 2.8 with a mean value of 628. Full CI re- 
sults were calculated using MOLPRO^iii^^'^ We see that 
the match between the MCCI and the FCI is both qual- 
itatively correct and quantitatively appears very good in 
sharp contrast to the incorrect behaviour of Hartree-Fock 
here (Fig. [15]). The good match of the MCCI with the 
FCI is shown in the NPE value of 0.628 kcal/mol despite 
the wavefunction only being a few-hundredths of a per- 
cent of the size of the FCI SD space. We note that in 
Ref. [33 the NPE, when using a slightly different basis set, 
was found to be 0.653 kcal/mol for ic-MRCCSD. 
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FIG. 15. MCCI {cmin = 10"=*), FCI and Hartree-Fock energy 
(Hartree) against the reaction coordinate x (Bohr) for BeH2 
using the cc-pVDZ basis. 

We see in Fig. [16] that the difference between the FCI 
and MCCI energy is smallest at x = A Bohr and increases 
as the reactants approach (x decreases) until about x — 
2 Bohr from then on it decreases as x decreases. We 
note that the maximum error does not correspond to the 
transition state. 



F. Ammonia inversion 



E. BeH2 

The model reaction for the formation of BeH2 was put 
forward in Ref. 33: to investigate CCSD. There CCSD was 
found to describe the system well even when the system 
had multireference character with two states. Recently, 
BeH2 has been considered as a test system for internally 
contracted multireference coupled cluster (ic-MRCC)24. 
In the aforementioned work the beryllium atom is at the 
origin and the hydrogen atoms have, in Bohr, the coordi- 
nates X, y and a;, —y where y = 2.54 — 0.46a; and x G [0, 4]. 
We run the calculations with no frozen orbitals and note 
that the full CI space is around 4 million SDs when ne- 



We investigate the ability of MCCI to reproduce the 
potential curve of NH3 as its trigonal pyramid inverts by 
moving through a planar structure. We use the cc-pVDZ 
basis and freeze one orbital. The NH bond length is 1.025 
angstrom and the hydrogens are at 120 degrees to each 
other. The NH bond makes an angle 9 with a line passing 
through the nitrogen and perpendicular to the hydrogen 
plane. 

The MCCI curve is displayed in Fig. [TTland compared 
with Hartree-Fock and FCI results which were calculated 
using MOLPRO.i^^^i^ The plot is mirrored around the 
line X = 90. We see that both MCCI with Cmin = 10"^ 
and Hartree-Fock recover the shape of the FCI curve with 
the minima and transition point in approximately the 
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FIG. 16. Error (Hartree) in the MCCI (cmi„ = 10"^) energy 
compared with FCI against the reaction coordinate x (Bohr) 
for BeH2 using the cc-pVDZ basis. The NPE is 0.63 kcal/mol. 

correct place. The MCCI curve appears to be a little 
too flat as it approaches the planar structure in that the 
energy does not change as much as the FCI or Hartree- 
Fock between 85 and 90 degrees. Despite this, for the 11 
points from 40 to 90 degrees, the NPE is 2.4 kcal/mol 
while Hartree-Fock has an NPE of 9.9 kcal/mol. The 
number of CSFs ranged from 1, 226 to 1, 824 with a mean 
value of 1, 629. This is in contrast to the size of the FCI 
SD space when neglecting symmetry of ^ 4 x 10*. 
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FIG. 17. MCCI {cmin = 10"^), FCI, and Hartree-Fock 
energy (Hartree) against angle (6) for NH3 using the cc-pVDZ 
basis and one frozen core orbital. 

The difference in energy between the MCCI and FCI 
calculation is displayed in Fig. [181 For this system we see 
that the error generally decreases as the planar transition 
state is approached. 

G. Hydrogen lattice 

We now consider a linear chain of hydrogen atoms. 
Such a system would perhaps be expected to be more 
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FIG. 18. Error (Hartree) in the MCCI (c^in = 10~^) energy 
when compared with the FCI against angle (6^) for NH3 using 
the cc-pVDZ basis and one frozen core orbital. The NPE is 
2.4 kcal/mol. 



suitable for modelling using techniques appropriate 
for strongly-correlated one-dimensional lattice systems 
such as density-matrix renormalization group (DMRG)~ 
methods which can in principle scale linearly with the 
size of the system; a linear 50 Hydrogen system has been 
considered using the DMRG for quantum chemistry in 
Ref. HI with a ST0-6G basis. The ID hydrogen chain 
therefore poses an interesting challenge for MCCI par- 
ticularly as the correlation progresses from dynamic to a 
large amount of, what could be considered, static corre- 
lation as the distances between the atoms increase: the 
FCI results using MOLPRCUiSS^ for a chain of 12 hy- 
drogens in the ST0-6G basis have only two coefficients 
greater than 0.05 ( -0.110 and 0.981) at 1.0 Bohr sepa- 
ration compared with 19 at 4.2 Bohr where the largest is 
only 0.22. 

In Fig. [Tni we see that for a chain of twelve hydrogens 
the Hartree-Fock energy increases much too fast as the 
distance between hydrogen atoms increases. MCCI, with 
Cmin = 10~^, is a fairly good match with the FCI re- 
sults. We note that there was a numerical issue with 
the MCCI results for R > 3.2 Bohr when using MOs. 
Here MCCI produced an energy slightly lower than the 
FCI energy. However when using orthogonal atomic or- 
bitals we found a MCCI energy of -5.694 Hartree at 4.2 
Bohr compared with the FCI result of -5.699 Hartree, 
this seemed to suggest that the use of very poor MOs in 
MCCI for this challenging system can lead to numerical 
problems when using CSFs in MCCI. The full configu- 
ration space consists of 853, 776 determinants while the 
MCCI wavefunction ranged from 342 at i? = 1 Bohr to a 
maximum of 6, 477 by i? = 2.8 Bohr. The mean value was 
2, 701 CSFs. The NPE for the displayed results was 15.8 
kcal/mol. In Fig. [20] we display the difference between 
the MCCI and FCI energies. We see that the MCCI en- 
ergy error increases from its value at i? = 1 Bohr until 
R = 2.8 Bohr. 
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FIG. 19. MCCI (cmin = 10"^), FCI, and Hartree-Fock 
energy (Hartree) against distance between atoms R (Bohr) 
for a chain of twelve hydrogen atoms with a STO-6G basis. 
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FIG. 20. Error (Hartree) in the MCCI (c^in = 10"^) energy 
compared with the FCI against distance between atoms R 
(Bohr) for a chain of twelve hydrogen atoms with a STO-6G 
basis. The NPE is 15.8 kcal/mol. 



We did not satisfactorily calculate a potential curve 
(Fig. I2ip for the 50 hydro gen system compared with 
DMRG based results of Ref. |33 which we suggest as due 
to the extremely large configuration space (10^®) of which 
we only included up to ~ 80,000 states (circa 10~^^% of 
the full space) which even assuming that the fraction of 
states needed for an accurate potential curve decreases 
with the size of the space still appears to be much too 
small a sample. We note that the DMRG results used lo- 
calized orbitals while the MCCI potential curve used HF 
MOs. We can see in Fig. [21] that the MCCI curve im- 
proves upon that of Hartree-Fock in this basis and retains 
the correct equilibrium bond length but also still displays 
incorrect behaviour at large R. This suggests that it does 
not sufficiently account for the static correlation for this 
system at the value of Cmin used (5 x 10~^). We note that 
the MCCI and Hartree-Fock curve for i? < 1.8 Bohr ap- 
pear to have moderately good agreement with the DMRG 



results which may be expected as the correlation is es- 
sentially dynamic here. 

As R becomes larger we would expect the system to be 
closer to a collection of non-interacting hydrogen atoms 
so the atomic orbitals might be expected to be a more 
efficient basis than the restricted Hartree-Fock molecu- 
lar orbitals. We create orthogonal atomic orbitals using 
the Gram-Schmidt procedure starting with the leftmost 
atomic orbital. The use of CSFs allow us to put one elec- 
tron in each atomic orbital and still have a correct S' = 
spin function. At R = 4.2 Bohr we find that we have 
an energy which is 0.4 Hartree above the DMRG result, 
but the MCCI program does not seem able to easily im- 
prove upon this single CSF as no other configurations are 
then found. Such an approach would not work well for 
slightly smaller R where both dynamic and static cor- 
relation is important — for R = 3.2 Bohr this approach 
yields a difference of 1.2 Hartree — and suggests that a 
larger fraction of the configuration space would need to 
be explored and that approximate natural orbitals could 
be useful. 
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FIG. 21. MCCI (cmin = 5 X 10"''), Hartree-Fock, 

LDMRG(500) (see Ref. Isj), and MCCI {c^,^ = 5 x lO-*) 
using orthogonal atomic orbitals (AOs). Energies (Hartree) 
against distance between atoms R (Bohr) for a chain of fifty 
hydrogen atoms when using a STO-6G basis. 



H. Ethylene torsional angle 

We now look at ethylene for which we use Rc-c — 

I. 325 angstrom, Rc-h = 1.090 angstrom and ZHCC — 
120.252 degrees. We use a cc-pVDZ basis with Cj^in = 
10^^ and vary the torsional angle. We find the tor- 
sional barrier in ethylene to be 75.52 kcal/mol which com- 
pares relatively favourably with the cc-pVDZ results from 
Ref. [H. There the barrier was found to be 68.2 kcal/mol 
when using GASSCF(12,12) while CASPT2 lowered this 
to 65.5 kcal/mol. The Hartree-Fock barrier was as high 
as 111.8 kcal/mol. We note that it is not clear whether 
exactly the same geometry was used and also whether it 
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was kept fixed, as in our results, or made more realistic by 
being allowed to relax at 90 degrees in the quoted litera- 
ture values (i.e., fully geometry optimizing the transition 
state structure). 

The number of CSFs increased from around 5, 900 to 
about 11,200 as the torsional angle increased to ninety 
degrees. The mean number of CSFs was 8, 250 which is 
8.3 X 10-12% of ^Yie size of the lO" SD FCI space when 
neglecting possible spatial symmetries. In Fig. [52] the re- 
sults have been mirrored about the line x = 90. The lack 
of a cusp at 90 degrees when using multireference per- 
turbation but its occurrence when using single reference 
methods such as CCSD was noted in Ref. js^ and we see 
here that a cusp appears to be avoided as the gradient of 
the MCCI curve decreases as it approaches 90 degrees. 
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FIG. 22. Energy against torsional angle (degrees) for ethylene 
with a cc-pVDZ basis and Cmm = 10"^. 



of the systems considered in this work. Quantifying the 
size-consistency by comparing the MCCI energy at the 
largest bond length and lowest cut-off used with the sum 
of the MCCI energy of the fragments reveals a difference 
of less than 10 kcal/mol for most of the systems con- 
sidered and less than 1 kcal/mol for those undergoing a 
single hydrogen dissociation, however N2 in a cc-pVTZ 
basis had a difference of 41.8 kcal/mol. 

We note that as the space became larger the fraction 
used tended to decrease, for example the results for the 
nitrogen and fluorine molecules in Table IIIIl which we 
would hope may be a general feature. We acknowledge 
that the NPE was larger for the nitrogen and fluorine 
molecules when the space became larger but compar- 
ing, e.g., nitrogen with BH we see that the NPE is only 
around 2.5 times larger but the fraction of the space is 
around 100 times smaller while the reduction in the space 
is even more pronounced for fluorine and the increase in 
NPE a little smaller. We note that the NII3 FCI has 
more than six times the states of BH but MCCI uses a 
fraction around 100 times smaller to give a slightly lower 
NPE. We also saw that nitrogen in a cc-pVDZ basis used 
2, 854 on average of a possible 10^ states while in a cc- 
pVTZ basis without frozen orbitals we could produce a 
well-behaved curve using just 10, 600 states, at the small- 
est cut-off considered, of the much larger 10^^ FCI space, 
although an NPE value was not available in the latter 
case. We acknowledge that a flxed value of Cmin restricts 
the size of the MCCI space so that the maximum fraction 
possible must decrease as the FCI space enlarges. How- 
ever our MCCI results are well within this limit and we 
note that the number of MCCI states appears to increase 
much more slowly than the FCI space. 



IV. SUMMARY 

We have seen that potential energy curves for small 
systems may be calculated by MCCI to approaching 
chemical accuracy, when compared with full configura- 
tion interaction results, even when the system is known 
to be multiconfigurational. Table IIIIl shows how NPE 
values of often around a few kcal/mol can be achieved 
using MCCI to construct a wavefunction that uses only 
a very small percentage of the size of the FCI space. For 
all the systems but BH, which has fairly similar NPE and 
mean single-point energy error, and the chain of twelve 
hydrogens, which had the largest fraction of states but 
worst NPE, we achieved a better NPE than the mean 
single point error demonstrating that MCCI can indeed 
be used for calculating a potential curve to higher accu- 
racy than its constituent single-point energies. We note 
that in general for a fixed Cmin > MCCI would not 
be expected to be size-consistent nor size-extensive, but 
the potential curves calculated here suggest, that for a 
low enough c,nin, MCCI is sufficiently size-consistent to 
correctly describe the approach to dissociation of many 



TABLE III. Table showing the mean CSFs to FCI SDs 
(without symmetry considerations) ratio, mean CSFs to FCI 
CSFs (without symmetry considerations), Non parallelity er- 
ror (NPE) (kcal/mol) and Mean single-point error (MSPE) 
(kcal/mol), at the smallest cut-off used, for the systems in- 
vestigated for which FCI results were available. 

System % FCI SD space % FCI CSF space NPE MSPE 



BH 


0.007% 


0.021% 


2.6 


2.1 


HE 


0.014% 


0.056% 


1.3 


2.6 


CH4 


0.012% 


0.049% 


0.6 


4.1 


C2 


0.003% 


0.013% 


4.9 


6.0 


F2 


8.3 X 10'^% 


4.9 X 10-*^% 


6.2 


9.0 


N2 


6.6 X 10"^% 


3.2 X 10-*% 


6.6 


11.9 


H12 


0.32% 


1.2% 


15.8 


6.0 


BeH2 


0.016 % 


0.052% 


0.63 


0.94 


NH3 


4 X 10"^% 


1.7 X 10"^% 


2.4 


8.4 



However, even with the possibility of a decrease in the 
fraction of the FCI space required as the system became 
very much larger, we found that a chain of fifty hydro- 
gens was to prove too challenging for MCCI. Here we saw 
that although the equilibrium geometry was correct this 
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was not an improvement on the Hartree-Fock results and 
for longer bond lengths, where static correlation would 
be considered important, the curve, although lower in 
energy than that of Hartree-Fock, increased too quickly. 
We suggest that this is due to the very large size of the 
configuration space (~ 10^^) and the very challenging 
system for standard quantum chemistry methods partly 
due to its strongly correlated nature: restricted Hartree- 
Fock is a very poor approximation at large bond lengths 
here and H12 required the largest fraction of the FCI 
space and gave the largest NPE. We note that as a sin- 
gle CSF using orthogonalised atomic orbitals was much 
closer to the DMRG energy here then the results could 
perhaps be improved by the use of approximate natural 
orbitals and allowing the algorithm to explore more of 
this space but this strongly correlated system may still 
require many states and seems to be more amenable to 



modelliug UHUig other lueUiodH. — Finally, Lhe poLeiiLial 
curve fui Lhu isumeiisaLiun uf ethylene was found Lu be 
relatively smooth and the calculated barrier compared 
fairly favourably with other computational results. Al- 
though we calculated potential surfaces, for which FCI 
results exist, to relatively high accuracy using just a very 
small fraction of the FCI space, we have noted two possi- 
ble limitations of MCCI when applied to potential energy 
surfaces. The first is that when the cut-off is not suffi- 
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